Definitions of entanglement entropy of spin systems in the valence-bond basis 
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The valence-bond structure of spin-1/2 Heisenberg antiferromagnets is closely related to quantum 
entanglement. We investigate measures of entanglement entropy based on transition graphs, which 
characterize state overlaps in the overcomplete valence-bond basis. The transition graphs can be 
generated using projector Monte Carlo simulations of ground states of specific hamiltonians or 
using importance-sampling of valence-bond configurations of amplitude- product states. We consider 
definitions of entanglement entropy based on the bonds or loops shared by two subsystems (bipartite 
entanglement). Results for the bond-based definition agrees with a previously studied definition 
using valence-bond wave functions (instead of the transition graphs, which involve two states). For 
the one dimensional Heisenberg chain, with uniform or random coupling constants, the prefactor of 
the logarithmic divergence with the size of the smaller subsystem agrees with exact results. For the 
ground state of the two-dimensional Heisenberg model (and also Neel-ordered amplitude-product 
states), there is a similar multiplicative violation of the area law. In contrast, the loop-based entropy 
obeys the area law in two dimensions, while still violating it in one dimension — both behaviors in 
accord with expectations for proper measures of entanglement entropy. 
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I. INTRODUCTION 

The concept of quantum mechanical entanglement is of 
broad interest in physics. One widely used quantitative 
measure of entanglement is the von Neumann entangle- 
ment entropy S'vn- Partitioning a many-body system in 
a pure quantum state into two contiguous subsystems, 
S'vn is defined as the von Neumann entropy of the re- 
duced density matrix of either one of the two subsystems. 
The insight that S'vn typically scales proportionally to 
the boundary area (the area lawi) was first developed 
in the context of black holes.— The area law has now 
become a key benchmark for characterizing states of in- 
terest in quantum information theory"^'** and condensed 
matter physics;^ An important application of this con- 
cept is to construct computationally tractable variational 
ansatze for ground states based on the area lawi^ 

While SvN is a clear and well established measure 
of entanglement, it is difficult to compute for strongly- 
correlated quantum systems in dimensions higher than 
one, where exact diagonalization and density matrix 
renormalization group approaches become inefficient j^ii^ 
Alternative definitions of entanglement entropy (or, more 
precisely, alternative measures of entanglement) are 
therefore also actively investigated. The Renyi entropies 
Sn (n = 1, 2, . . .) are often usedfii and have the appeal- 
ing property that Si = Svn- Recently it was realized^^ 
that S2 (and in principle also S„ for n > 2) can be com- 
puted for quantum spin systems using projector quan- 
tum Monte Carlo (QMC) simulations in the valence bond 
(VB) basis. ^- Previously, a measure Svb of entanglement 
entropy explicitly making use of the VB basis was also 
introduced^'*'^''' within the context of this QMC method. 
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Generalizing an exact result for a single VB stati 

a superposition, Svb is given by the average number 



of VBs connecting two subsystems. While defined ex- 
plicitly using the VB basis, this quantity also can be 
evaluated using the density matrix renormalization group 
(DMRG) method^^ and is, therefore, not completely tied 
to the VB basis. Its asymptotic behavior has also been 
found exactly using analytical methods for the Heisen- 
berg chaiuj^ 

In this paper, we formulate a different measure Syg of 
entanglement entropy for quantum spin systems in terms 
of the transition graphs characterizing overlaps of VB 
basis states. The transition graphs are generated in pro- 
jector Monte Carlo simulations of the ground state of 
a hamiltoniau) ^'^d^ or in Monte Carlo sampling of bond 
configurations of variational states such as the amplitude 
product statesJ^i^ Like the previous definition of VB 
entanglement entropy, which we henceforth call Sys i the 
transition-graph definition SyB involves VBs shared by 
two subsystems, but the weighting is different because 
the bond configurations in the projected bra and ket 
states are sampled using their individual wave functions 
and overlap (which depends on the number of loops in 
the transition graph). We show that Syg scales with 
the subsystem size in accord with an exact result^^ for 
the previous definition Syg in the one dimensional (ID) 
Heisenberg chain. The corrections to this form are much 
smaller than in the previous definition, however. Thus, 
the different weighting procedure appears to reduce the 
subleading size corrections. 

Both Syg and Syg violate the area law in the case 
of the Neel-ordered ground state of the two-dimensional 
(2D) Heisenberg model. We here argue that this is be- 
cause definitions based on single VBs typically will over- 
estimate the entanglement, due to the over-completeness 
of the VB basis. To remedy this, we propose an alter- 
native measure Sioop of entanglement entropy based on 



the loops of the transition graphs, i.e., S'loop is the aver- 
age number of transition-graph loops shared by the two 
subsystems. Loops correspond to maximally entangled 
groups of spins, and the number of loops shared by the 
subsystems is therefore an appropriate measure of entan- 
glement. We show that S'loop of the 2D Heisenberg model 
(and also in a generic variational amplitude-product state 
with Neel order) obeys the area law and has an additive 
logarithmic correction, in contrast to the multiplicative 
logarithmic corrections affecting the bond-based defini- 
tions. Thus, it appears that the loop entropy scales in the 
same way as the Renyi entropy (as observed in Ref. 1121 ) 
and the standard von Neumann entanglement entropy 
(where one would expect such scaling) and, thus, may be 
a convenient (more easily computable) stand-in for these 
definitions. 

The outline of the rest of the paper is as follows. 
In Sec. [n] we summarize the properties of the VB ba- 
sis needed for our definitions and calculations, and also 
briefly review amplitude-product states;^ the VB pro- 
jector QMC method^^ii^ and the loop-gas picture^ii^ 
of VB states. The VB and loop entropy definitions are 
discussed and tested in Sees. IIIII and IIVI We conclude in 
Sec. |V] with a summary and discussion. 



II. VALENCE-BOND BASIS AND METHODS 



connects sites on the two sublattices) with all positive 
coefficients;^- 
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With the singlet state of two spins i,j denoted by (i, j), 
a bipartite VB basis state is defined as 
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where sites i and j are on different sublattices and each 
site appears exactly once in the product. Thus, there 
are {N/2)\ different basis states \v) which form an over- 
complete basis in the singlet subspace. 

The overlap {v'\v) between two VB basis states can 
be expressed in terms of transition-graph loops, as illus- 
trated in Fig. [TJ Each lattice site is connected to one 
bond from \v) and one from (w'|, and all the bonds there- 
fore form closed loops. The matrix element {v'\v) is a 
product of factors arising from these loops. To evaluate 
{v'\v) it is more convenient to work in the standard spin-z 
basis and rewrite a VB state \v) as 
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We will consider the spin-1/2 Heisenberg hamiltonian, 
written in the form 



H — — 2_^ JijSij, {Jij > 0), 



(1) 



where {i,j) denotes nearest-neighbor spins on a lattice 
with periodic boundaries and Sij is a singlet projector. 



S,y = 1/4 - S, ■ Sj 



(2) 



and Jij > are antiferromagnetic coupling constants. 
We will study ID chains with uniform and random cou- 
plings, as well as uniform 2D square lattices. 

Here, in Sec. Ill Al we discuss the VB basis in which we 
carry out all calculations. In Sec. Ill B] we discuss varia- 
tional VB amplitude-product states, which give some im- 
portant insights into various types of ground states and 
their entanglement properties. The projector QMC cal- 
culations that we use for unbiased calculations are briefly 
reviewed in Sec. Ill CI The loop-gas picture of VB states 
was flrst introduced by Sutherland i^iSS In Sec. IIIDI we 
discuss it in a somewhat broader sense, which we will rely 
on for the definition of loop entropy in Sec. IIVI 



A. The valence-bond basis 

The ground state of i? on a bipartite lattice with 
an even number N of spins is a total-spin singlet and 
can be expanded in bipartite VB states (i.e., each bond 



where a labels the spin states that are compatible with 
the valence bonds in \v) (i.e., fj, or 4,t spin configurations 
on each bond) and ng-f- is the number of f spins on sub- 
lattice B, the sign following if the singlet ^ is defined 
such that i and j are on sublattices A and B, respec- 
tively (which corresponds to Marshall's sign rule^^ for a 
bipartite system) . From the orthonormality of the ordi- 
nary basis of states {IS*^, Sf, • • • , S^j)}, the only nonzero 
terms in {v'\v) expressed using (O are those with spin 
configurations common to both the VB states \v) and 
\v') . This corresponds to spins forming staggered pat- 
terns, tit • • • i or tti ■ • • tj around each loop. There are 
two such staggered configurations of each loop and the 
signs in ([5]) cancel in the overlap, which, thus, is given 
by2i 
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where A^o is the number of loops in the transition graph. 
This is illustrated in Fig. [TJ 

In the VB description many physically quantities are 
determined by the statistical properties of the transition- 
graph loops. For instance, the normalized matrix element 
needed for computing the spin-spin correlation function 
is given by 
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where [i,j] and [i][j] denote sites i and j belonging to the 
same loop and different loops, respectively, and the sign 
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FIG. 1: (Color online) Transition graph (right) determining 
the overlap of two VB states \v} and \v'}. Note that in the 
restricted VB basis considered here, all bonds connect sites on 
different sublattices (with sites on sublattices A and B here 
indicated by darker and lighter circles) . The cluster of spins 
defined by a loop can be in two states, with all spins t on 
sublattice A and J, on B, or vice versa. 



in the case [i,j] is + and — for spins on the same and 
different sublattices, respectively. 

Note that the transition graph loops define clusters of 
spins that do not necessarily resemble geometric loops on 
the lattice, because the VBs can be of any length and the 
bonds forming a loops can cross each other (as seen in 
Fig. H]). Note also that the staggered spin configuration 
along a loop also corresponds to staggering of all spins in 
the loop in the sense of the two sublattices, i.e., all spins 
on sublattice A are parallel and opposite to those on B. 
This is the origin of the signs in ([7]). 



B. Amplitude-product states 



the amplitudes can be variationally optimized. iSiSS, Cal- 
culations for the Neel-ordered ground state of the 2D 
Heisenberg model suggest that the fully optimized am- 
plitudes decay as a power-law, /i(|rij|) = l/|r,y|P, with 
p « 3.^^'^^ There is reason to believe that the decay ex- 
ponent in fact is exactly p = 3, as this is the exponent ob- 
tained in an analytical mean- field-like treatment. ^^ Some 
aspects of the critical ID chain can also be captured with 
amplitude-product states<^ 

The properties of quantum systems are often modi- 
fied dramatically by introducing quenched (static) ran- 
domness: e.g., quantum phase transitions with disorder 
can lead to new universality classes. For the Heisenberg 
chain, it is known that any amount of quenched random- 
ness will drive the system into a state well approximated 
by a random singlet state — a single VB basis state (which 
is of the "nested" type, with no crossing VBs) with arbi- 
trary VB lengths obtained according to a strong-disorder 
renormalization-group (SDRG) scheme^rJfl In this case, 
the transition-graph loops coincide exactly with the VBs, 
i.e.. No = N/2 in ([6]), and many asymptotic properties 
follow directly from the length distribution of the VBs. 

Amplitude-product states (as well as the special case 
of the SDRG random singlet states) are useful variational 
states and we will consider some aspects of their entan- 
glement properties. In many cases completely unbiased 
results are needed, however. One way to achieve this is 
with projector QMC simulations in the VB basis, which 
we briefly discuss next. 



In an amplitude-product state, ^'^ the expansion coef- 
ficients in ^ are products of amplitudes h{rij) > (to 
maintain Marshall's sign rule for the ground state of a 
bipartite model). 
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with the vector r^ denoting the "shape" (the lengths in 
all lattice dimensions) of bond {i,j) in \v). This form 
applies to a translationally invariant system, while in a 
non-uniform system one should use h{ri,rj) depending 
on both end-points of the bonds. 

For a given set of amplitudes, the expectation value of 
some operator O can be written as 

,^X ^ (V^0|O|^0) _ E.„t,' Wyy- "^^^,1^" 

where the configuration weight is given by 

Ww = A„A„- {v'\v) = A„A„-2^""^/^ (10) 

The ratio {v' \0\v) / {v' \v) can normally be related to the 
loop structure of the transition graph,— i^ e.g., Eq. ([7]) 
in the case of a spin correlation function. 

The expectation value (O is ideally suited for evalu- 
ation using Monte Carlo sampling methods /^'^" and all 



C. Projector QMC method 

In the VB QMC methodi^i^ the ground state of the 
hamiltonian ([1]) is projected out stochastically, by apply- 
ing a high power of H to some trial singlet state Itpt) in 
the VB basis; [~HY^ |i/'t) —J' IV'o) for a large m (up to an 
irrelevant normalization). A good trial state, such as an 
optimized amplitude-product state for a uniform system 
or a single VB state obtained with the SDRG procedure 
for a ID random chain, can be used to optimize the con- 
vergence properties of such a scheme, but the final result 
is not sensitive to \^t) as long as m is sufficiently large 
(i.e., the method is unbiased). 

There are two formulations of the VB QMC method^^ 
generating either the ground state wave function or the 
ket and bra versions of the ground state needed to eval- 
uate expectation values. In the former case, stochastic 
application of the projector ff"* (for sufficiently large m) 
on the trial state produces VB basis states distributed 
proportionally to the expansion coefficients A,j in ([3|), 
i.e., these coefficients are not known (and would be much 
more complicated than the simple amplitude products) 
but importance-sampling of the contributions to ff™|'0i), 
which resemble terms in a path integral, generate them 
probabilistically (with many paths contributing to a sin- 
gle coefficient A^,). 

In the "double projection" method, an expectation 



value is formally given by Eq. ([9]), but the expansion 
coefficients A^, are again not known. The sampling of 
paths obtained from the combined ket iJ™ | ■(/;() and bra 
(-0t| ^™ states, i.e., (V't| -ff^™ IV't): leads to a series of VB- 
pair configurations distributed according to the weights 
VF™- in ([T0|. 

In both the wave-function and expectation-value pro- 
jection schemes, one can employ efficient loop updates for 
generating the states and transition graphs, as discussed 
in detail in Ref. [l^ In the sampling procedures of such 
algorithms, one uses the combined spin-bond basis and 
applies a high power if™ or H^™ in a way similar to the 
" operator- loop" update in the finite temperature stochas- 
tic series expansion QMC methodi^ The computational 
effort scales as 0{m), and to converge calculated quanti- 
ties to the ground state to has to be scaled as ex A^°, with 
a typically in the range 1 — 2 (depending on the model's 
low-energy energy spectrum and the quality of the trial 
state IV't))'^ 



D. Loop-gas description 

The loop-gas view of VB sampling was suggested by 
Sutherland some time ago<2L2^ Consider the two VB con- 
figurations in a transition graph {v'\v) with the associ- 
ated weight ([TU)) for an amplitude-product state. In that 
case, the factor XyXy' does not depend on details of the 
bond arrangements, only on the total number iVf,(r) of 
bonds of all shapes r in the two states \v) and \v') (i.e., 
no bond correlations are included). In addition to the 
factor 2 in the overlap ([6|) for each loop in the transition 
graph, for each loop containing more than two bonds we 
can consider swapping all bonds belonging to a given loop 
between \v) and \v'), as illustrated in Fig. [21 This swap- 
ping does not affect the weight (|10p of the joint configu- 
ration of the two states (although the weights Ay and X^i 
of the individual states are affected) . The two bond con- 
figurations and two staggered spin patterns for each loop 
corresponds to a loop fugacity 4 for all loops of length 
4 or larger. The shortest, length-2, loops have fugacity 
2, as the bond swapping in this case does not affect the 
bond configurations. 

Instead of sampling bonds weighted according to (TTOl) , 
one can think of sampling loops with the weight 



W = 2^°^^h^°^>^^Y[h{r] 



Nb(r) 



(11) 



where iVo(2) and No{> 2) denote the number of loops of 
length 2 and larger than 2, respectively, and the unimpor- 
tant factor 2~^/^ in ^ has been omitted. The product 
of amplitudes can here be thought of as originating from 
the shapes of the loops. In special cases, such as An- 
derson's resonating valence-bond (RVB) state including 
only the shortest bonds (of length r = 1),—"— the weight 
only depends on the number of loops. This was the case 
considered by Sutherland, ^^ who also studied generaliza- 
tions of the loop gas in which the loop fugacities, Z2 = 2 





FIG. 2: (Color online) Two transition graphs {v'\v) with the 
same loops structure (with four loops, labeled 1,2,3,4). The 
two cases shown here illustrate how bonds within a loop (here 
loop 4) can be swapped between \v) and {v'\ (represented by 
solid and dashed bonds, respectively). This kind of bond 
reconfiguration leaves the weight unchanged in the case of an 
amplitude-product state, but not in a more general state in 
which there are correlations between the bonds. 



and Z>2 = 4 in (jlip . can take arbitrary values and, thus, 
phase transitions can be studied as a function of these 
fugacities. 

The loop-gas weight function pT|) does not apply to 
states beyond the amplitude-product description. In the 
exact ground state of a given hamiltonian, one would in 
general expect bond correlations. The product of two 
VB expansion coefficients in (jlOp then changes when an 
intra-loop bond swap of the type illustrated in Fig. [2] is 
carried out, thus invalidating the form (jlip . One can 
still, however, sum up the weights obtained as a result of 
all such bond swaps which leave the loop structure intact, 
and this way, in principle (but hardly in practice), write 
down a weight for a loop configurations which is more 
complicated than (ITU) . The loop configurations gener- 
ated in the QMC projector method represent a stochas- 
tic implementation of this more general loop-gas picture. 
Note that even in this generalization, each loop is asso- 
ciated with a factor of two arising from the two allowed 
staggered spin configurations of the cluster of sites de- 
fined by the loop. 



III. VALENCE-BOND ENTANGLEMENT 
ENTROPY 



A singlet state formed by two S = 1/2 spins is maxi- 
mally entangled and has the maximum value of the von 
Neumann bipartite entanglement entropy; 5'vn = 1 (mea- 
sured in bits). The simplest case for computing 5vn of 
a many-body system is a single VB basis state \v). For 
a given bipartition, its entanglement entropy is just the 
number of singlets hab connecting subsystems A and 
B.^^ It is not easy to compute 5vn foi" a superposition 
of VB states, however. In Refs. 14 and 15 a straight- 
forward generalization of the result for a single VB state 
was proposed as an alternative definition of entanglement 
entropy for an arbitrary superposition of VB states, using 
the average number of bonds connecting the two sub- 
systems. With the VB projector method for the wave 
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FIG. 3: (Color online) (a) VB entropies for subsystems of 
size ^ in a Heisenberg chain of length L = 2^^ as a func- 
tion of the logarithm of the conformal distance; log2(^') = 
log2[(L/7r) sin(7r£/L)]. The lines have slope 7 = 41n(2)/7r^, 
as in the exact asymptotic form.— (b) Disorder-averaged VB 
entropies for the ID random chain, as a function of log2(/). 
The line has slope 7 = ln(2)/3. 



function, one stochastically evaluates 






^VB 



(12) 



Here we investigate a different definition, using the ex- 
pectation value 



'S'vB = (»^Ab), 



(13) 



evaluated using Eq. ^ with O — tias. ■ Here it should be 
noted that we define tiab as just counting of the bonds 
crossing the two subsystems, and that this counting can 
be performed in the transition graph of {v'\v) either in 
\v) or {v'\. For a given configuration these numbers are 
different (i.e., the operator is not hermitian), but the 
averages are the same. In Sec. |V] we will address fur- 
ther potential problems in interpreting (tiab) as a stan- 
dard quantum mechanical expectation value. For now, 
we consider this quantity as an interesting aspect of the 
transition graphs, the statistical properties of which are 
uniquely determined for a given hamiltonian with short- 
range interactions (as we will further argue in SecFV]). 

For a single VB state S'vb — 5'vb' but in a superpo- 
sition these quantities are different, because of the over- 
completeness and associated different weighting of the 
VB configurations. We next present projector QMC re- 
sults for both 5'yB and S'yg in ID and 2D Heisenberg 
systems and discuss their scaling as a function of the size 
of the smaller subsystem of the bipartition. 



A. Homogeneous chain 

We consider first the standard Heisenberg chain with 
uniform interactions. For a large segment of size £, em- 



bedded in a chain of length L, the von Neumann entan- 
glement entropy has the asymptotic behavior 5'vn(^) = 
(c/3) log2 {(.') + So, where ^ = [L/tt) sin(7r^/L) is the con- 
formal length, c = 1 is the central charge, and sq is a 
non-universal constanti^ The VB entropy S'yg is known 
to diverge in the same way, but with a different factor, 
7 = 41n(2)/7r^ w 0.281.^^ Previous calculations are con- 
sistent with 7 < 1/3 for large chains fiS but, as can be 
seen in Fig. |3l[a) for a chain of 2^^ spins, there are still 
large non-asymptotic corrections which make it difficult 
to verify the factor preciselyjii^, On the other hand, our 
results for Syg are completely consistent with the known 
7 over a large range of subsystems. This may appear 
surprising, because the exact calculation is based on the 
wave- function definition S'yg ,— not the expectation value 
Syg. It is plausible, however, that Syg and S'yg scale 
asymptotically with the same 7 (with different additive 
constants), but that S'yg is less affected by subleading 
scaling corrections. 



B. Disordered chain 

Now we turn to the disordered chain, with random 
couplings generated from the uniform distribution in the 
interval (0, 1]. The ground state of this system is known 
to be well approximated by a random singlet, where 
all spins form a single non-crossing VB state with ar- 
bitrary bipartite bond lengthsi^^ Using an SDRG anal- 
ysis, Refael and MooreiS showed that the disorder aver- 
age of the von Neumann entanglement entropy SvN in 
such a state scales logarithmically with a universal co- 
efficient 7 = ln(2)/3. This result should hold also for 
the VB entanglement entropy of the Heisenberg chain 
when L — > 00, although the ground state is not exactly a 
single VB basis state — there are fluctuations around the 
dominant SDRG configuration.'^^ Our results for both 
S'yg and S'yg, shown in Fig. [3l^b), are consistent with 
7 = ln(2)/3 (considering some remaining finite-size ef- 
fects for the largest subsystems). 



C. Two-dimensional system 

We next consider the 2D Heisenberg model, which we 
have studied on L x L lattices with L up to 256. Fig. 0] 
shows results for S'vb(^) of square subsystems of linear 
size i. The data converge very rapidly with Lior £ < L/4. 
There is clearly a multiplicative logarithmic correction 
to the area law, Syu{£)/£ ex log2(^), as found previously 
for 5yB in Refs. [l0lll4l Both VB entropy definitions, 
thus, violate the area law in this case. The von Neumann 
entropy, on the other hand, should obey the area law, 
S'vn(^) oc £, although, because of difficulties in calculating 
this quantity in an unbiased way, it has not been possible 
to confirm it unambiguously.!^ The recent calculation of 
the Renyi entropy ^2 is, however, in agreement with the 
area lawJ^ 
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FIG. 5: (Color online) The loop entropy S'loop of a subsystem 
(gray square) is defined as the average number of loops in the 
transition graph shared by the subsystem and the rest of the 
system. In this example two loops are shared. 



FIG. 4: (Color online) The VB entropy Svb divided by the 
linear dimension ^ of a square sub-block in 2D L x L lattices 
with L = 32 - 256. 



We conclude that Syg has an advantage over 5'yg, in 
the sense that its logarithmic prefactor in ID converges 
faster to the result of the exact calculation in Ref. [l^ (al- 
though we do not know whether the faster convergence 
is generic). Neither of these definitions serves well as a 
proxy for the von Neumann entanglement entropy 5vn, 
however, because the area law is violated in 2D by a 
multiplicative logarithmic correction. Only additive cor- 
rections to the area law are expectedi^i^ 



IV. LOOP ENTANGLEMENT ENTROPY 



We will here consider a measure of entanglement en- 
tropy iSioop in terms of shared loops in the transition 
graph, as illustrated in Fig.[SJ Since the number of shared 
loops must be smaller or equal to the number of shared 
bonds, we have S'loop < •S'yg, and it is therefore clear 
that S'loop cannot be a better approximation to S'vn than 
Syg in the case of the ID Heisenberg chain. We will be 
mainly interested in S'loop of the 2D system, but will also 
investigate it in ID. 

Before discussing the actual definition of Sioop further, 
let us consider for a moment the wave function. Instead 
of a bond-singlet pairing of the spins, one might regard a 
state as a superposition of products of loop-cluster states, 



w„ 



i*) = E(8)i^«)- 



(14) 



The loop-gas picture discussed in Sec. IIIDI suggests 
a potential reason why the VB entanglement entropies 
(IT^ and ((T5)) violate the area law in 2D systems: The 
VB basis consists of singlet pairs, but when consider- 
ing an expectation value, constraints related to the over- 
completeness are imposed on the spin states beyond the 
singlet pairing. These constraints are similar to multi- 
spin entanglement, as there are two spin states for each 
loop (the two staggered spin configurations on the loops) 
and these loops are, thus, analogous to maximally entan- 
gled sets of spins (albeit in an overlap matrix element, 
not a wave function). All bonds therefore do not carry 
a full unit of entanglement entropy (because they are 
entangled with other bonds in the same loop), and the 
VB entropy S'yg may therefore overestimate the actual 
entanglement entropy. This could also be the case with 
Syg, because it, too, is based on a superposition of states 
with only pair-wise entanglement (and, as we showed in 
Sec. mil Syg and Syg have the same scaling proper- 
ties). On the other hand, in ID Syg and Syg actually 
underestimate the entanglement entropy (in relation to 
the von Neumann entanglement entropy, which is slightly 
larger than the bond entropies in this case). The intu- 
itive picture of entanglement entropy in terms of bonds is 
therefore, as a consequence of the overcompleteness, not 
always quantitatively correct. 



where Na is the number of loops in component a and an 
individual cluster state of n spins is of the form 



IV'(O). 
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(lti;2---;n>L-i;it2---uL). (i5) 



This is a maximally entangled state of two staggered spin 
configurations along the loop of n spins (or, equivalently, 
all even numbered spins are on sublattice A and all odd 
ones on i3, or vice versa). With such as loop-cluster wave 
function, it is natural to regard a loop shared by two sub- 
systems in a bipartition as carrying one unit of entangle- 
ment entropy. Thus, we define the looj) entanglement 
entropy as 



Sloop — (A as), 



(16) 



where Aab counts the number of loops shared by sub- 
systems A and B (i.e., loops passing through both the 
subsystems) in a given bipartition. 

The loop-cluster view of entanglement entropy is re- 
alized in SDRG calculations for the random transverse- 
field Ising model)^ where the clusters have parallel spins 
and the ground state is given by a single product of such 
cluster states. The generalization to a superposition of 
cluster products has the same motivation as the general- 
ization of the entanglement entropy of the single VB state 
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FIG. 6: (Color online) Loop entropy of Heisenberg chains 
with (a) uniform and (b) random couplings (versus the log- 
arithms of the conformal and linear subsystem sizes, respec- 
tively). Both lines have slope 7 = ln(2)/3. 



appearing in the SDRG scheme for the random Heisen- 
berg chaini^ to a superposition. l^iiS, In general, it is not 
known, however, how to write the wave function of a 
Heisenberg system as a superposition of cluster states 
(and clearly there is no unique way of doing so, as such a 
basis is massively overcomplete). As we have seen, such 
loop-cluster superpositions do appear in the loop-gas pic- 
ture, and it is then natural to define Sioop as above using 
the transition graph loops, as illustrated in Fig. [SJ We 
will explore this measure of entanglement here. 

One may again question whether an expectation value 
such as (Aab) is a bona fide quantum mechanical expec- 
tation value. We will discuss this further in Sec. |V] and 
here only consider it, like Sy^, as a statistical property 
of the transposition graphs generated in the projector 
QMC scheme (or in sampling of an amplitude-product 
state, which we will also investigate). 

Clearly, the loop entropy is a boundary property, and 
we can write either S'ioop(^) or Sioop{B) for a biparti- 
tion {A,B). It is easy to demonstrate the sub-additive 
property 5'ioop(^i U A2) < 5'ioop(^i) + 5'ioop(^2) for two 
different bipartitions, {Ai,Bi) and (^2,-62)- The essen- 
tial properties of an entropy are thus satisfied. For a 
single VB configuration, loops and bonds coincide and, 
thus, S'loop — Sy-g ~ Syg = S'vN in this case. 



A. ID and 2D Heisenberg models 

We first consider the loop entropy for ground states 
of Heisenberg models obtained by the projector QMC 
method. Figs. [6] and [7] show results for ID and 2D sys- 
tems, respectively. Interestingly, in ID, the behavior is 
consistent with a logarithmic divergence with a prefactor 
7 — ln(2)/3 « 0.231 (to within a statistical error of about 
1%) for both the uniform and random chains (with uni- 
formly distributed couplings between and 1), i.e., the 
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FIG. 7: (Color online) The loop entropy divided by the linear 
dimension ^ of a square sub-block in several L x L lattices. 
The curve shows a fit of the L = 256 data {£ > 2^) to the 
form Sloop = jJ.£ + alog2{£) + so, with /i ~ 0.51, a ~ 1.12, and 
So ~ —0.34. The quality of this fit is statistically sound, with 
X^/dof < 1 when ^ > 8 data are included. The error bars are 
much smaller than the symbols, of the order 10~* for Sioop/£- 



same as the bond-based entropy for the random chain 
(Fig. [3|). The result is expected for the random chain, 
because the asymptotic SDRG state of such a system is 
a single VB state, in which S'loop = S'vb = '5'ioop {^^d this 
is the exact result), but it is curious that the same value 
obtains for the uniform chain as well. This may also be 
taken as a flaw of the loop-based entropy (if one wants a a 
quantity mimicking the von Neumann entropy as closely 
as possible), because it is further from S'vN (7 = 1/3) 
than the bond-based estimates (7 = 41n(2)/7r^ sa 0.281). 
Nevertheless, it is encouraging that the logarithmic di- 
vergence is still captured. 

In the 2D system, S\oo-p/ (■ converges to a finite value 
with increasing (, and L. Note the rapid convergence 
as a function of the full system size in Fig. [T) For the 
largest system (L = 256) the results are described very 
well by an area law with an additive logarithmic correc- 
tion; S'loop = M^ + afog2(^) + Sf). The area law prefactor 
/i ss 0.51 is, intriguingly, in good agreement with Svn/^ 
obtained by Kallin et alt^ for wide ladder systems [ac- 
counting for the ln(2) included there in the definition 
of Svn]- In this case the loop entropy is, thus, a better 
stand-in for Svn (which is expected to obey the area law) 
than the bond-based quantities. 

Note that the presence of long-range Neel order in 
2D implies the existence of system-spanning loops in 
the transition graph. ^^ This explains why Sioop is much 
smaller than Sye in Fig. [5] (i.e., because of the large aver- 
age loop size, the number of loops in the transition graph 
is much smaller than the number of bonds). 



B. Neel-ordered 2D amplitude- product state 

We also investigate the scaling behavior of the loop en- 
tropy for a 2D amplitude product state, ^ with the ex- 
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FIG. 8: (Color online). Loop entropy of the amplitude- 
product state with amplitudes decaying with the bond length 
r as h{r) = 1/r^, which is the asymptotic form in the opti- 
mal amplitude product state for the 2D Heisenberg model. ^^ 
The entropy divided by the system length saturates for large 
systems, indicating applicability of the area law (with an ad- 
ditive logarithmic correction described by the same form as 
in Fig. [T] a fit is indicated by the solid curve). 



pansion coefficients given by ([S]) with h{r) — l/r^. This 
corresponds to the asymptotic form of the optimal varia- 
tional amplitudes for the 2D Heisenberg model suggested 
by previous calculations J^i^ We here do not use fully 
optimized amplitudes (which show deviations from the 
1/r^ for short bonds), because our aim is to investigate a 
generic Neel-ordered state, to compare with the results of 
the specific case of the ground state if the 2D Heisenberg 
model in Fig. [T] 

The sublattice magnetization in the amplitude-product 
state with h{r) — 1/r^ for all r is m<; « 0.27 (extracted 



from the system size dependence of 



for systems of 



size L up to 256), somewhat below the known value rus ~ 
0.307 for the Heisenberg model. The results in Fig. [5] 
show a behavior very similar to the Heisenberg ground 
state result in Fig. [71 with a correction to the area law 
which can be described as an additive logarithm. Because 
of the lower value of the sublattice magnetization, the 
average loop size is smaller, and, thus, the number of 
loops in the system (including those shared by the two 
subsystems) is larger, leading to higher overall value of 
S'loop than in Fig. [5l 



V. SUMMARY AND DISCUSSION 

Using the transition graphs characterizing state over- 
laps in the valence bond basis of S* = 1/2 spins, we 
have explored two measures of bipartite entanglement 
entropy, S'yg and S'loop- The former extends the defini- 
tion Syg of valence-bond entropy (Refs. IlilllSl ) based on 
shared valance bonds in the wave function to the tran- 
sition graph, while the latter is based on shared loops 
(motivated by the loop-gas picture^iiS^ of spin systems). 
Using an efficient loop algorithmj^^, we were able to ob- 



tain unbiased QMC results for these quantities in large 
systems; 2^'^ and 2^^ spins in ID and 2D, respectively. 

In the Heisenberg chain, Sy^ exhibits a logarithmic 
divergence, with a pref actor which agrees well with the 
exact value'^^ 7 = 41n(2)/7r^ already for small subsys- 
tems. In contrast, observing this scaling with the wave- 
function definition Sy^^^' requires very large systems, 
due to significant subleading correctionsiiii 

For the Neel state of the 2D Heisenberg model, both 
VB entropies violate the area law, exhibiting multi- 
plicative logarithmic corrections J^ We have argued that 
single-bond definitions typically overestimate the amount 
of entanglement, because of the over-completeness of the 
VB basis. The loop definition S'loop exhibits only an ad- 
ditive size correction to the area law in 2D, in agreement 
with general expectations for standard definitions J^i^i^ 

Important relationships have been established in recent 
years between the subleading behavior of the entangle- 
ment entropy, topological order, and quantum-criticality. 
For instance, the subleading term in a 2D gapped sys- 
tem is a constant determined by the quantum dimen- 
sion of the excitations of the topological phase. "^'^ For 
critical systems in the universality class of z = 2 con- 
formal quantum-critical points, there is a universal ad- 
ditive logarithmic subleading term, which depends only 
on the shape of the subsystem partition and the central 
charge:^ An additive correction to the area law in the 2D 
Neel state was also found for the Renyi entropy S2, but 
the systems were too small to determine its asymptotic 
form, I- We have shown that the loop entropy Sioops cap- 
tures the essential desired features of an entanglement 
entropy, scaling logarithmically with the subsystem size 
in ID and obeying the area law for the 2D Neel state. 
Being relatively easy to calculate, Sioops offers opportu- 
nities to study various aspects of entanglement entropy 
on large spin lattices in other situations of great interest, 
e.g., at unconventional quantum-critical points4ii^ 

In Sec. lIIII we already commented on the fact that prop- 
erties such as Syg and Sioop that are defined using spe- 
cific geometrical properties of the transition graphs (i.e., 
not following directly from a given operator acting on the 
spins) in the valence bond basis are not necessarily well 
defined expectation values of some hermitian operators. 
Indeed, it has recently been shown that the entanglement 
entropies defined in this way are dependent on exactly 
how a state is represented in the overcomplete valence 
bond basis4^ This may suggest that these quantities are 
ill-defined. However, given a hamiltonian H, the projec- 
tor QMC method generates the transition graphs in a 
unique way, independently of the trial state used (which 
we have also verified explicitly). 

The valence-bond projector technique itself is closely 
tied to the completely generic "loop-operator" represen- 
tation of the path integral for a singlet state^^i^ and 
therefore the statistical properties of the transition graph 
loops (including Syg and Sioop) are not really tied to the 
particular QMC scheme, only to the valence-bond basis. 
The definitions are tied to the hamiltonian, in the sense 



that H generates the transition graphs. We can there- 
fore not, in general, evaluate S'yg and i'loop uniquely 
just based on an arbitrary state, but first need to find 
the "parent hamiltonian" of the state (and furthermore, 
that parent hamiltonian should have short-range interac- 
tions only, so that it is unique — support for this generic 
statement is discussed in Ref. 1451 ). and use it to gener- 
ate the transition graphs. States defined based on "un- 
biased" bond superpositions (which includes contribu- 
tions from all different ways of expressing the state in 
terms of singlets obeying Marshall's sign rule), such as 
the amplitude-product states (and perhaps generaliza- 
tions of them including bond correlations) can also be 
studied, as we have done here in a 2D case (and where 
it is important to note that for such a state with Neel 
order, we obtained results in good agreement with the 
projected ground state of the 2D Heisenberg model). In 
practice, we expect Syg and S\oo-p to be useful primarily 
in QMC studies of ground states of specific hamiltonians. 
In the loop-operator formulationfi^i^ one can think of 
clusters of spins (defined by transition-graph loops) as 



forming dynamically in imaginary time. The entangle- 
ment entropies S'yg and 5'ioop correspond to entropy mea- 
surements averaged over equal-time "snapshots" of en- 
tangled clusters in this time evolution. Exactly how this 
dynamic aspect relates to the standard wave-function 
picture of entanglement entropy is not clear at present, 
but the results obtained in this paper suggest that there 
should be a close relationship between them. 
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